Proton— deuteron elastic scattering above the deuteron breakup 
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The complex Kohn variational principle and the (correlated) hyperspherical harmonics method 
are applied to study the proton-deuteron elastic scattering at energies above the deuteron breakup 
threshold. Results for the elastic cross section and various elastic polarization observables have 
been obtained by taking into account the long-range effect of the Coulomb interaction and using 
a realistic nucleon-nucleon interaction model. Detailed comparison to the accurate and abundant 
elastic proton-deuteron experimental data can now be performed. 



o\ : 

A number of nucleon-nucleon (NN) potentials are now 
available which can be used in a non-relativistic approach 
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to understand nuclear structure. Those usually referred 
to as realistic potentials accurately fit the data base se- 
lected by the Nijmegen group Q with a x 2 per datum 
close to one. They produce quite reasonable values for 
the binding energies of different light nuclei with small 
deviations from the experimental values. One way of 
reducing these deviations is to add three-nucleon (3N) 
interaction terms determined in a semiphenomenological 
j — ■ way. 

When such potentials are used to calculate scattering 
states, some observables, such as the elastic N — d dif- 
ferential cross sections, are well predicted (2|^]. How- 
ever, the situation is different for a number of observ- 
ables which strongly depend on the nuclear interaction 
in specific waves. Examples of these are spin-dependent 
observables, such as the vector or the tensor analyzing 
1 1 ■ powers. For these quantities there are significant dif- 
ferences between theoretical estimates and experimental 
values. This is a strong signal that there are deficiencies 
in the theoretical models adopted. Detailed theoretical 
and experimental investigations are therefore necessary 
Jv>( , in order to understand the reasons for this problem. 

In this respect, the study of the p — d scattering pro- 
cess is of particular relevance since very accurate mea- 
surements exist for a large set of observables and kine- 
matical regimes. The Faddeev theory has been very suc- 
cessfully applied to the n — d process (|,|||| but the 
extension to the p — d case, taking properly into ac- 
count the long range Coulomb repulsion, presents a num- 
ber of difficulties which have been the object of exten- 
sive researches ^,0]. Accurate calculations of TV — d 
scattering below the deuteron breakup threshold (DBT) 
have also been performed in the frame of the so-called 
pair-correlated hyperspherical harmonic (PHH) expan- 
sion technique ||. The incorporation of the Coulomb 
potential in the PHH approach does not present particu- 
lar difficulties, and p — d scattering observables have been 
calculated using realistic NN+3N potentials (see refs. ||] 
and (^]). The technique used is variational and based 
on the use of the Kohn variational principle (KVP). The 
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extension to the study of the four-body (p — 3 He and 
n — 3 H) zero energy scattering has been given in ref. 

In ref. jlO| the authors used the complex form of the 
KVP to describe n—d and p—d scattering above the DBT, 
using a semi-phenomenological nuclear s-wave potential, 
while the Coulomb interaction was included without any 
partial wave projection. For the n — d case, the results 
obtained were in close agreement with the benchmarks 
obtained solving the Faddeev equations in configuration 
and momentum space 

In the present paper the study of N — d scattering 
above the DBT is extended to the case of realistic NN 
potentials. Cross sections as well as vector and tensor 
polarization observables for n — d and p — d scattering 
for nucleon incident energies up to 10 MeV have been 
calculated using the NN Argonne AV18 potential |fL2| . 
These results are compared to the available experimental 
data @,|l5| . Some of the n — d results are compared with 
the values obtained by the Bochum-Cracow group ||. 

The details of the variational approach used by us can 
be found in refs. §Q|^. The applicability of the KVP 
above the DBT when the Coulomb interaction is taken 
into account deserves some attention. Here a brief discus- 
sion of its validity for the description of elastic scattering 
is given. A more general discussion will be reported else- 
where ]Ti| . The scattering wave function (w.f.) W is 
written as a sum of two terms. The first term, ^>c, de- 
scribes the system when the three-nucleons are close to 
each other. For large interparticle separations and ener- 
gies below the DBT, it goes to zero, whereas for higher 
energies it must reproduce a three outgoing particle state. 
It is written as a sum of three Faddeev-like amplitudes 
corresponding to even permutations of the particle in- 
dices 1, 2, 3. Each amplitude 4'c'(x,, yj), where Xj,y^ 
are the Jacobi coordinates corresponding to the i-th per- 
mutation, has total angular momentum JJ Z and total 
isospin TT Z and it is decomposed into N c channels using 
the LS coupling, namely 

*c(xi,yi) = y^<j> a (x h yj)y a (jk,i) (1) 

a. 

y a (jk, i) = { [Y la (xi)Y La (fe)] Aa [si k sl] Sa }^ [t{%] TTz , (2) 

where Xi, yi are the moduli of the Jacobi coordinates and 
y a is the angular-spin-isospin function for each channel. 
The two-dimensional amplitude <p a is expanded in terms 
of the PHH basis 



<t>a{xi,yi)=P la+Lm ~ 5/2 f*(zi) 



^u<* K { P ) {2) Pt L ° 



L K 



(3) 



where the hyperspherical variables are defined by the re- 
lations Xi — p cos <fii and yi = p sin <f>i, f a ( x i) is a P a h 
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correlation function and ^ P^ L ((f>) is a hyperspherical 
polynomial. 

The second term in the variational scattering w.f. de- 
scribes the asymptotic motion of a deuteron relative to 
the third nucleon. It can also be written as a sum of 
three amplitudes in terms of the ingoing and outgoing 
solutions of the asymptotic N-d Schroedinger equation. 

o+ Si/ ( Xi , yi ) = n%j(^, yi ) - £ J s s L s L :ni%,j(^, yi ) , 

L'S> 

(4) 

where J Sf i f t , are the elastic S'-matrix elements. 

The three-nucleon scattering w.f. for an incident state 
with relative angular momentum L, spin S and total an- 
gular momentum J is 

*LSJ= + n£ 5J (x i , y< )] , (5) 

i=l,3 

and its complex conjugate is ^2sj- ^ variational esti- 
mate of the trial parameters in the w.f. VP \ s j can be ob- 
tained by requiring, in accordance to the complex KVP, 
that the functional 

['sif,] = j s s L s L : + i(n S j\H - e \*l>s'j) > (6) 

be stationary. 

The validity of the KVP above the DBT and with 
charged particles for the elastic 5*-matrix elements is 
briefly discussed below. Let us consider the w.f. Wjjgj 
describing the p — d process for an energy E, and a trial 
approximation of it W^cr. Both w.f. can be written 
in the form given in eq. (|^), with the assumption that 
for the exact one the sum in eq. (Jfy is not truncated 
at any level. The hyperradial functions and S'-matrix 
coefficients entering the w.f. ^lsj will be specified by 
an overline to distinguish them from the corresponding 
trial quantities. In the asymptotic region p — > 00, the 
hyperradial functions are superpositions of ingoing and 
outgoing waves 

p l " +L "U K {p) - J^^^ 2 ^™' 13 ^^ 
a'K' 

_ J2(e- ixlos2Qp )™'A K ',e iQ <' , (7) 

a'K' 

where Q 2 = M^E/h 2 and the ^-matrix originates from 
the Coulomb potential. Since we are interested in the 
process p + d — > (p + d) + (p + p + n), the boundary 
conditions to be imposed are 

B% = 0, for all K, a . (8) 

For p — > 00 we can specify four regions, characterized 
by different ranges of values of the hyperangular variables 
[CI = 4>, x, y]. The [fi&] region is the breakup region where 
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all the particles are well separated. The region where the 
particles j and k are close each other, while particle i is 
very far from them, is hereafter denoted by [f2j]. There 
are three such regions, corresponding to the cases i = 
1, 2, 3. Let us consider the integral 

I = (*Isj\(H - m+ SIJ ) R - ^ usl j\(H - E)* + LS j)r 

(9) 

where stands for the integration in the six- 

dimensional volume with p < R (and R — > oo). Only 
the differential operators present in H contribute to i\ 
After integrating by parts, the contributions come from 
the hypersurface at p = R where the trial and exact wave 
functions have reached their asymptotic behavior. 

Let us write / = /;,+ Y]j—i Ii: where 1^ (Zj) is the 
contribution coming from the region [f2&] In 
the asymptotic functions are vanishingly small and 

lb reduces to 

I^Y^UKiP^^iP)-^)^^)] • (10) 



dp 



The above form has been obtained after orthonormalizing 
the PHH basis elements at p = oo. Using the asymptotic 
behavior given in eqs. (0,||) for both the exact and trial 
hyperradial functions, lb — > as R — » oo. In the three 
regions [Qi], the breakup part of the w.f. can be neglected 
since it gives contributions which go to zero as i? _3//2 , 

therefore Y^=i^i ^ J Sll< ~ J Sll'- Finally, using the 

fact that (H—E)i& LS j = 0, it is possible to show that the 

f s s' 

functional [ J >S£f/] differs from J S LL , only quadratically 

in the difference e = 9 — 

The crucial points of the proof are i) the outgoing 
boundary conditions satisfied by and ii) the null con- 
tribution to Ii of the breakup part. The presence of the 
Coulomb potential introduces a distortion in the outgo- 
ing waves which, essentially, does not change the main 
points of the demonstration of the KVP given for the 
elastic part of the S-matrix in the n-d case |l8| . 

The variation of the functional with respect to the hy- 
perradial functions leads to the following set of coupled 
equations 



E [ A $ 0°) -T2 + B w (p) -r + <%* (p)- 



^-EN^'(p)]u^(p) = D^ k (p) . (11) 

For each asymptotic state ^ 2S+1 ^Lj two different inho- 
mogeneous terms can be constructed corresponding to 
the asymptotic Q^gj functions with A = in, out. The 
numerical technique used to solve the above set of equa- 
tions imposing outgoing boundary conditions at a finite 



value of the hyperradius p = po is given in ref. |10|. Es- 



sentialy, the solutions of eq. (Mb for p > p are obtained 
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as series in lip imposing the outgoing boundary condi- 
tions of eqs.(f?],||). In the case of n — d scattering such 
solutions evolve as outgoing Hankel functions H^\Qp). 
In the region p < po the hyperradial functions have been 
expanded as 



M 
m=0 

(12) 

where z = "/p and 7 is a nonlinear parameter. The func- 
tions Lm (z) are Laguerre polynomials. The parameters 
A m and 7 arc determined by the variational procedure. 
The functions defined above are matched to the outgoing 
solutions at po. The value of the matching radius po is 
not critical and a value of po ~ 100 fm has been found 
to be satisfactory. 

The functions u a ,fc(/?) are the solutions of eq. (|l] 
where all the couplings between the differential equa- 
tions have been neglected (and applying outgoing bound- 
ary conditions). Their inclusion is necessary since the 
functions (p) show an oscillatory behavior already for 
p > 30 fm. To reproduce such a behavior would require 
a rather large value for M in eq. (|l2|). However, the in- 
clusion of the terms u a} K allows values of M similar to 
those needed for describing N — d scattering below the 
DBT @. 

In order to check the convergence properties of the 
PHH expansion for energies above the DBT we first 



solved the same problem treated in ref. 10 using the 
present technique. All phase-shift and inelasticity pa- 
rameters were reproduced with the same previous accu- 
racy, i.e. with a precision of four figures. As compared 
with ref. jl0| , the dimensions of the matrices involved in 
the eigenvalue problem came out reduced by one order 
of magnitude. 

Let us start studying N — d scattering above the DBT 
using the AV18 interaction. The nuclear elastic S-matrix 
has been calculated up to total angular momentum states 
J = 11/2+. This includes all partial waves with relative 
angular momentum L < 4. Higher partial waves (up to 
L = 8) were included in the calculation of the observables 
using the Born approximation ||. For each J T state all 
channels with £ a + L a < Kq have been included. The 
number of hyperradial functions has been increased until 
convergence was reached. The maximum value Kq = 6 
was found appropriate to obtain the elastic scattering 
observables within an accuracy of 1%. The pattern of 
convergence in terms of Kq was studied in ref. ]13| ] for 
energies below the DBT and a similar behavior has been 
observed here. 

High quality measurements of p — d scattering have 
been presented in ref. p^| . Cross sections and proton 
analyzing powers have been measured up to Ei a f, = 18 
MeV, deuteron analyzing powers and tensor analyzing 
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powers up to E [ab = 9 MeV [E d = 18 MeV). In fig. 1 
our theoretical predictions for these observables are com- 
pared to the data at Ei a \, = 5 MeV. In fig. 2 the same 
set of observables at P; a h = 10 MeV are compared to the 
data of ref. pi] . In addition to the p—d calculations (solid 
line), the n — d results (dashed line) are also shown for the 
sake of comparison. A good agreement between theory 
and experiment is observed for the differential cross sec- 
tion. The already known puzzle has been found again for 
the vector analyzing powers A y and iT\\ which are un- 
derpredicted by about 30%. The tensor analyzing powers 
are rather well described, with small underpredictions at 
the second minimum in T20, the second maximum in T21 
and the minimum in T22. These differences increase with 
energy. The origin of these discrepancies can be analyzed 
in terms of phase shift and mixing parameters. For ex- 
ample, in ref. phase-shift analyses were performed at 
Ei a b = 2.5 MeV and 3.0 MeV with the conclusion that 
small differences in the P-wave phase-shifts and mixing 
parameters were responsible for the discrepancy in the 
A y and iTu observables. This problem with the P-wave 
parameters seems to persist also at higher energies. The 
small discrepancies in the tensor observables could origi- 
nate from higher partial waves. In fact, the tensor observ- 
ables are particularly sensitive to phase-shift and mixing 
parameters for L > 2. At Ei a b < 3.0 MeV, just below 
the DBT, these parameters are small due to centrifugal 
barrier effects, but at the energies considered here their 
contribution becomes appreciable. 

Faddeev calculations in momentum space for n — d 
elastic scattering at P; a fc = 5 and 10 MeV have been 
presented in ref. Q for several potential models includ- 
ing the AV18 potential. Our corresponding results are in 
complete agreement with those reference calculations. 

In conclusion, p — d elastic cross sections and polar- 
ization observables have been calculated with a realistic 
interaction for energies above the DBT up to Ei a b = 10 
MeV and taking into account Coulomb interaction ef- 
fects. Accurate calculations of p — d observables and 
their comparison with the available experimental data 
may give stringent tests of the existing models of NN 
and 3N interactions. The extension of the present tech- 
nique to higher energies and to the breakup cross sections 
will be the subject of a forthcoming paper. 
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FIG. 1. Differential cross section da/dfl, proton analyzing 
power A y , deuteron analyzing power iTn and tensor analyz- 
ing powers T20, T21 and T22 calculated at Ei a b = 5 MeV and 
compared with the data of ref. jl4| (circles with error bars). 
The solid (dashed) lines are the p — d (n — d) results. 

FIG. 2. As in fig. 1, but for E lah = 10 MeV. The data 
reported here are from ref. ]l5(]. 
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